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Abstract 

The English grain aphid, Sitobion avenae (F.) (Homoptera: Aphididae), is a dominant and de- 
structive pest in wheat, Triticum estivum L. (Poales: Poaceae), production regions in China and 
other grain-growing areas worldwide. Patterns of gene expression of the S. avenae-resistant syn- 
thetic wheat line 98-10-35, the S. ave/iae-susceptible linel376, and their hybrid population, and 
the differences in segments between 98-10-35/1376 F3 resistant plants and resistant parents of 98- 
10-35, as well as those between the F3 resistant and susceptible populations, were examined with 
differential display reverse transcription PCR. The results showed that five patterns of differential 
expression were detected between the progeny and its resistant parents: 1) The gene was silenced 
in one of the parents; 2) Special expression showed in the progeny; 3) Expression was consistent 
with the resistant parents; 4) Up expression showed in the progeny but not in the parents; 5) 
Down expression showed in the progeny but not in the parents. Paired Mest results were not sig- 
nificant; however, the probability value (0.9158) indicated that gene expression on the RNA level 
were consistent with resistant bands found in F3 resistant individuals and resistant parents, as well 
as the F3 resistant and susceptible populations. For both the F3 of 98-10-35/1376 and the parents, 
the total number of amplified bands was 202, with an average of 25.3 per primer. The number of 
differential bands was 116, with an average of 14.5 per primer amplified and a polymorphism 
ratio of 56.3%. In the present study, differential expression genes in the resistant line 98-10-35 
were all up-regulated. Among them, gene expression of resistant groups in the F3 population was 
in agreement with patterns 2, 3, and 4. However, the susceptible line 1376 did not have this gene 
expression on the RNA level. This pattern is expected to be used to select and analyze target 
genes from the same F3 population and the resistant parents. The results suggest that it can be 
employed as a new method for molecular assisted breeding. 
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Introduction 



The English grain aphid, Sitobion avenae (F.) 
(Homoptera: Aphididae), is an important de- 
structive pest in wheat, Triticum estivum L. 
(Poales: Poaceae), production regions in Chi- 
na and the world (Alkhedir et al. 2010; 
Razmjou et al. 2011; Liu et al. 2012). In Chi- 
na, the influence of S. avenae is spread over 
an area of about 16.7 million hm in wheat 
production areas, especially in the Yellow 
River and Huai River basins, the North China 
Plain, and the southwest, northwest, and mid- 
dle Yangtze River basins in China (Zhang et 
al. 2009). As a result, S. avenae infestation 
has caused reductions in wheat grain yield by 
as much as 42% (George and Gair 1979; Xu et 
al. 1998). Chemical control has been used to 
try to solve this problem, but it has led to a 
severe environmental pollution (Flickinger et 
al. 1991; Daily et al. 1998). However, effi- 
cient, large-scale control measures are not 
available at present. Alternatively, it has been 
shown that plant resistance is the most com- 
mon approach in integrated pest management 
programs to counter the effects of S. avenae 
(Razmjou et al. 201 1). 

In recent years, the effect of wheat variety on 
major population parameters of S. avenae- 
susceptible or S. avenae-resistmit wheat varie- 
ties have been extensively investigated by 
numerous scientists worldwide (Ozder 2002; 
Delp et al. 2009; Dogimont et al. 2010). Hu et 



al. (2004) studied German and Chinese wheat 
cultivars to determine their physical resistance 
locus to S. avenae by using an electrical pene- 
tration graph technique. Previous studies 
typically focused on the mechanism of differ- 
ential expression levels of inducible resistance 
to wheat germplasm (Zhu et al. 2005; 
Dedryver et al. 2010). Some differential ex- 
pression genes were obtained from these 
studies, but the results could not fully explain 
the mechanism of resistance of wheat grain to 
aphids (Zhu et al. 2005; Zhao et al. 2009; 
Zhang et al. 2012). Despite the progress in 
this area, the exact molecular mechanisms of 
the differential expressions of wheat genetic 
resistance to S. avenae require further study. 
Differential expression of constitutive re- 
sistance of wheat varieties to S. avenae is still 
unknown. 

It has been suggested that breeding and culti- 
vating aphid-resistant germplasm is the most 
economical, safe, and effective way to man- 
age S. avenae infestation. Breeding plants that 
are resistant to S. avenae is considered an im- 
portant method for the integrated pest 
management. In recent years, many effective 
resistance genes have been isolated with the 
rapid development of biotechnology, which 
enables molecular-assisted breeding of re- 
sistant wheat to be used in integrated pest 
management. 
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Differential display reverse transcription PCR 
(DDRT-PCR) has been widely used in studies 
of signal transduction and morphological de- 
velopment of plants and in identifying and 
cloning plant resistance genes associated with 
stress and disease. Many plant genes and 
cDNA fragments have been obtained using 
this technique (Xiong et al. 2002; Chen et al. 
2006; Delp et al. 2009). This technique shows 
potential for studying the functional genomics 
in plants. However, no study on the differen- 
tial expression resistance-related genes in 
wheat parents, the F3 generation, or the back- 
cross populations has been reported. In 
addition, the expression patterns between the 
F3 resistant populations and parents or sensi- 
tive plants have not been examined. 

The wheat line 98-10-35 showed resistance to 
S. avenae (Du et al. 1999; Hu et al. 2011; 
Wang et al. 2011). The resistance was con- 
trolled by a single dominant gene (Wang et al. 
2011). The S. avenae-YQsistmit line 98-10-35, 
the S. avenae-svLSceptible line 1376, the F3 
generation, and the backcross population of 
98-10-35/1376 were investigated in this study. 
The objective of this study was to examine the 
expression pattern of the F3 resistant popula- 
tions, the differences in expression between 
parents and F3 resistant populations, and the 
differences between the sensitive and resistant 
plants. Expression patterns could provide an 
identification method for selecting target 
genes from candidate genes. 

Materials and Methods 

Wheat lines 

The wheat lines used in this study were 98-10- 
35, 1376, Amigo, and resistant and susceptible 
plants 98-10-35/1376 F 3 (177, 178) BCi and 
BC2 population. All the wheat materials were 
provided by Henan University of Science and 
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Technology, Northwest A & F University, 
China. 

Field trials 

Materials were planted on the experimental 
farm of Henan University of Science and 
Technology. The field trials were investigated 
over two successive seasons (2010/2011 and 
2011/2012). Each wheat germplasm line was 
sowed on two 100 cm long rows, 24 cm apart. 
Two replications were established with a 
completely randomized design. Guard rows 
were set, and susceptible line 1376 was plant- 
ed on them. No insecticide was applied during 
the entire growing season. 

Identification of resistance to aphids 

Resistance identification was conducted under 
the natural infection condition. For each 
germplasm line, the number of S. avenae on 
the 10 most seriously injured stems was 
counted to investigate the total number of S. 
avenae (four replicates for each germplasm) 
during the jointing stage and grain- filling 
stage. Infested plants were examined every 7 
days. The wheat grain aphid index (average 
number of a certain germplasm of S. avenae 
per plant/average number all germplasm of S. 
avenae per plant) was calculated. The wheat 
grain aphid index was summarized in seven 
scales, as shown in Table 1 (Painter et al. 
1958). 

Primers of DDRT-PCR 

Primers in this study were synthesized by 
Sangon Biotech Co. ( www.life-biotech.com ). 



Table I . The index of what germplasm resistance to Sitobion 
avenae. 



Resistance level 


Resistance 


Aphid index 


0 


Immunity 


0 


1 


Highly resistant 


0.01-0.30 


2 


Moderately resistant 


0.31-0.60 


3 


Lowly resistant 


0.61-090 


4 


Lowly susceptible 


0.91-1.20 


5 


Moderately susceptible 


1.21-1.50 


6 


Highly susceptible 


> 1.50 
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Table 2. The primers used in DDRT-PCR analysis (Lao 

2009). 


Name of primers 


Sequence base of primer 


H-T11A 


AAGCTTTTTTTTTTTA 


H-T11G 


AAGCTTTTTTTTTTTG 


H-T11C 


AAGCTTTTTTTTTTTC 


H-AP-1 


AAGCTTGATTGCC 


H-AP-2 


AAGCTTCGACTG T 


H-AP-3 


AAGCTTTGGTCAG 


H-AP-4 


AAGCTTCTCAACG 


H-AP-5 


AAGCTTAGTAGG C 


H-AP-6 


AAGCTTGCACCAT 


H-AP-7 


AAGCTTAACGAGG 


H-AP-8 


AAGCTTTTACCGC 






Figure I . The quality testing of total RNA (A) and the first- 
strand cDNA (B). 1-6 are resistant and susceptible parent 
plants and resistant and susceptible plants of F3 population and 
BCi population, respectively. M represents the DNA marker. 
High quality figures are available online. 



The primers consisted of three anchor primers 
and eight random primers, as shown in Table 
2 (Lao 2009). Meanwhile, the house-keeping 
gene 18srRNA primers were used to examine 
the quality of total RNA and successfully re- 
versed transcription. The sequences were 5'- 
CGTCCCTGCCCTTTGTACAC-3' and 5'- 
CTCCATGTC ATCCC AGTTG-3 ' . 

Synthesis of the first strand of cDNA 

For RNA extraction, 1 cm leaf samples were 
immediately submerged in 2 mL Trizol (Invi- 
trogen, www.lifetechnologies.com ), 
homogenized for 1 min with a tissue homoge- 
nizer, and stored at -80° C until further 
process (Tiangen Biotech Co. 
www.tiangen.com) . The total RNAs were re- 
versely transcribed to form the first chain of 
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cDNA by PrimeScript RT Reagent Kit 
(TaKara Biotechnology Co., 

www.clontech.com/takara ). Quality testing of 
RNA extraction and first strand cDNA are 
shown in Figure 1 . 

DDRT-PCR 

The 15 |iL reaction mixtures consisted of 1.5 
|iL reaction products of the first chain of 
cDNA, 1.5 |iL Mg-free lOxPCR buffer, 1.5 
unit of Taq DNA polymerase, 1.3 |iL 2.5 mM 
of MgCl 2 , 1.3 |iL 2.5 mM of dNTP, and 1.0 
|iL 10 |iM of each primer pair and ddF^O. Af- 
ter 5 min of denaturation at 94° C, 
amplifications were conducted for 25 cycles, 
with each cycle consisting of denaturing 30 
sec at 94° C, annealing 1 min at 45° C, and 
extension 1 min at 72° C followed by 15 cy- 
cles, each one consisting of 30 sec at 94° C, 1 
min at 50° C, 1 min at 72° C, followed by ex- 
tension at 72° C for 10 min. 

Data Analysis 

A matching Mest was done with SPSS 13.0 
(IBM, www.ibm.com ). The amplified bands 
were analyzed with Quantity One software 
(Bio-Rad, www.bio-rad.com ). The stable 
bands were recorded in two repeat amplifica- 
tions from 50 bp to 1000 bp. The lane with a 
detected band was marked T, otherwise it 
was marked '0'. 

Results 

Identification for resistance to aphid 

The results showed that wheat line 1376 was a 
highly susceptible line, and 98-10-35 and 
Amigo were moderately resistant lines. The 



Table 3. Resistance evaluation of parent or population to 


Sitobion avenae. 








No. 


Parent or population 


Number of plants tested 


Aphid index 


Resistance scales 


1 


98-10-35 


16 


0.465 t 0.092 


MR 


2 


Amigo 


16 


0.4047 : 0.167 


MR 


3 


1376 


16 


3.504 ± 0.056 


HS 


4 


98-10-35/1376 BC, 


16 


0.481 ±0.096 


MR 


98-10-35/1376 BC, 


1.519 ±0.102 


HS 


5 


98-1 0-35/1 376F 3 


71 


0.428 ± 0.074 


MR 


98-10-35/ 1376 F, 


2.014 ±0.275 


HS 


6 


178 


16 


1.0520 ±0.0823 


LS 


7 


177 


16 


0.4670 ±0.1495 


MR 
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Figure 2. 8% non-denatured PAGE picture of differential 
results. I — 1 0 and M represent two F3 plants of Amigo/I 376, 
1 77 line, 1 78 line, F 3 plant of 98- 1 0-35/ 1 376, 1 376, Amigo, 98- 
IO-35, BCi plant of 98- 1 035/ 1 376, BC 2 plant of 98- 1 0- 
35/ 1 376, and marker, respectively. High quality figures are 
available online. 




Figure 3. Expression patterns of hybrids and their models in 
wheat. High quality figures are available online. 



P F: 



Figure 4. Gene expression models of parent and its F3. Each 
type of three bands from left to right represents male parent 
1 376, F3 plant of 98- 1 0-35/ 1 376, and female parent 98- 1 0-35, 
respectively. High quality figures are available online. 



BCi and F3 populations were segregated in 
two parts, including moderately resistant and 
highly susceptible. Average aphid indexes of 
resistant and susceptible plants are listed in 
Table 3. 
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Analysis of the differential transcript- 
derived fragments 

Three anchor primers and eight random pri- 
mers consisted of 24 DDRT-PCR primer 
pairs, which were applied to test all plants in 
order to find the differential expression tran- 
script-derived fragments (TDF) (Figure 2). 

Differences of expression patterns in re- 
vealing a candidate gene 

Twenty-four pairs of primers were employed 
to conduct the repeating test to screen wheat 
resistance genes to S. avenae in all samples. 
The results showed that five patterns of differ- 
ential expression were detected between the 
progeny and its resistant parents (Figure 3): 
(1) The gene was silenced in one of the par- 
ents (a); (2) One of the parents appeared in 
bands, while others had no bands. Special ex- 
pression showed in progeny (b); (3) The gene 
only occurred in hybrid offspring bands, and 
two parents had no bands. Express consisten- 
cy showed with resistant parents (c), which 
means that bands for resistant parents and off- 
spring appeared, but did not appear in the 
sensitive parents; (4) Up-expression only 
showed in progeny and not in either parent 

(d) ; (5) Both the parents and offspring ex- 
pressed the resistance bands, but the amount 
of expression in the offspring was much more 
than the parents. Down expression only 
showed in progeny and not in either parent 

(e) . 

According to the amplification, five kinds of 
differential expression patterns can be detect- 
ed in the samples. Analysis and comparison of 
the amplified bands showed that there were 
slight differences among the F3 generation, 
backcross offspring of 98-10-35/1376, and the 
parents. In the genetic process, its expression 
model mainly appeared in five patterns, as 
shown in Figure 4. 
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Table 4. Results of TDFs of candidate genes. 





Primer wquencc (5'-3') 


No. ofampUned 
cONA 


Differential cDNA 


No. ofF3 
ampllfled cDNA 


No. of 98-10-35 
polymorphic cDNA 


No. of 1376 
polymorphic cDNA 


Number 


% 


H-TIIA 


aagcii i i n ri itia 


30 


16 


53J 


12 


9 


9 


H-AP-I 


AAGCTKiATTGCC" 


H-TIIA 


AAGC I 1 1 1 I 1 1 1 1 I IA 


32 


14 


43.8 


11 


9 


0 


H-AP-2 


AA(X"TTC"(iAl"TtiT 


H-TIIA 
H-AP-4 


AAOCTTCK AA( (i 


30 


21 


70 


13 


IS 


6 


H-T11Q 




30 


21 


70 


13 


8 


2 


H-AP-2 


AAUCTTCUACTOT 


H-riui 




20 


7 


35 


3 


2 


1 


H-AP-4 


AAGCTTCTCAACU 


H-TI1C 




25 


13 


52 


6 


9 


0 


H-AP-2 


A AGt TI C( i AC TG I 


H-TIIC 




10 


S 


SO 


5 


0 


H-AP-< 


AACKTTAGTAGG C 


HTIIC 


AAGC I I I 1 I I I I n k: 


25 


19 


76 


11 


9 


11 


H-AP-7 


AAGCTTAACCiAGG 


Total 


202 


116 


450.1 


74 


62 


29 


Mean 


25.3 


M J 


56 1 


9.3 


7.8 


3.6 



/ \ 

Table 5. Results of differential fragments of candidate genes. 



Primer 


Primer sequence (5'-3') 


No. of amplified 
bands 


Differential cDNA 


Number 


% 


H-TIIA 


AAGCTTTTTTTTTTTA 


28 


17 


60.7 


H-AP-1 


AAGCTTGATTGCC 


H-TIIA 


AAGCTTTTTTTTTTTA 


27 


21 


77.8 


H-AP-4 


AAGCTTCTCAACG 


H-T11G 


AAGCTTTTTTTTTTTG 


23 


10 


43.5 


H-AP-2 


AAGCTTCGACTGT 


H-TIIA 




28 


21 


75 


H-AP-2 


AAGCTTCGACTGT 


H-T11G 


AAGCTTTTTTTTTTTG 


23 


9 


39.1 


H-AP-4 


AAGCTTCTCAACG 


H-T11C 


AAGCTTTTTTTTTTTC 


20 


10 


50 


H-AP-2 


AAGCTTCGACTGT 


H-T11C 


AAGCTTTTTTTTTTTC 


6 


4 


66.7 


H-AP-5 


AAGCTTAGTAGG C 


H-T11C 


AAGCTTTTTTTTTTTC 


32 


19 


59.4 


H-AP-7 


AAGCTTAACGAGG 


Total 


187 


111 


472.2 


Mean 


23.4 


13.9 


59 





Expression differences analysis of TDFs in 
F 3 of 98-10-35/1376 and parents 

TDF groups for F 3 of 98-10-35/1376 and par- 
ents were amplified via eight pairs of 
differential DDRT-PCR primer. Analysis and 
comparison of amplification figures showed 
some differences in TDFs among the F3 of 98- 
10-35/1376, backcross offspring, and parents. 

TDFs changed to a certain extent in type and 
number on F3 of 98-10-35/1376 and parents, 
with the statistics segment size being 50-1000 
bp. Statistical results of differentially ex- 
pressed fragments are shown in Table 4. 

The results showed that the total number of 
amplified bands was 202 for the F3 of 98-10- 
35/1376 and parents, with an average of 25.3 
bands per primer pair. There were 116 differ- 
ential bands, with an average of 14.5 bands 
per primer pair and a polymorphism ratio of 
56.3%. From the analysis and comparison of 
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mutual bands on F3 of 98-10-35/1376 and par- 
ents, there were 74 differentially expressed 
bands for F3 generation of 98-10-35/1376, 
with an average of 9.3 bands per primer. 
There were 62 differentially expressed bands 
for the S. avenae-resistant line 98-10-35, with 
an average of 7.8 bands per primer. There 
were 29 differentially expressed bands in 
1376, with an average of 3.6 bands per primer. 
The fragments matching Mest showed that p = 
0.9158 and showed no difference for the F3 
generation of 98-10-35/1376 and 98-10-35 
differentially expressed. The F3 of 98-10- 
35/1376, the line 98-10-35, and line 1376 ex- 
hibited significant differences in the amplified 
products, i.e., gene expression products only 
existed in the F3 generation. 

Different TDF displays of F 3 generation re- 
sistant and sensitive single plants 

TDF groups for resistant single plants and 
sensitive single plants of F3of 98-10-35/1376 
were amplified by using eight pairs of differ- 
ential primers with a TDF segment size of 50- 
1000 bp. 

The types and numbers of TDFs of the re- 
sistant single plants and sensitive single plants 
varied greatly on the F 3 of 98-10-35/1376. The 
results showed that the total number of bands 
amplified was 187 in the resistant single plant 
group and sensitive single plant of F3 of 98- 
10-35/1376, with an average of 23.4 bands per 
primer (Table 5). There were 111 differential 
bands, with an average of 13.9 bands per pri- 
mer amplification and a polymorphism ratio 
of 59%. The differential expression gene 
products appeared in the F3 of resistant and 
sensitive single plants. The resistant single 
plants and sensitive single plants of the F3 
population of 98-10-35/1376 showed differ- 
ences in amplified products when different 
primers were used. 
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In summary, differential expression genes in 
resistant line 98-10-35 were all up-expressed. 
The gene expression of resistant groups in the 
F3 generation was in agreement with the pat- 
terns 2, 3, and 4 (Figure 2). The sensitive 
cultivar 1376 failed to express. 

Discussion 

DDRT-PCR 

At present, there are several methods for iden- 
tifying differentially expressed sequence tags, 
such as northern blotting, differential screen- 
ing, microarray analysis, subtractive 
hybridization, serial analysis of gene expres- 
sion, cDNA amplified fragment length 
polymorphism, and DDRT-PCR (Lao 2009). 
This study focused on the molecular mecha- 
nism of wheat germplasm resistance to S. 
avenae. From the aspect of gene expression, 
differential expression was explored by breed- 
ing resistant cultivars to aphid, and some of 
available genes were obtained by using 
DDRT-PCR. 

DDRT-PCR is a cDNA differential display 
technique that combines RNA reverse tran- 
scription with traditional PCR technology. 
This technique is highly sensitive and has a 
wide range of applications, such as testing the 
differential expression of genes at the cellular 
level (Xiong et al. 2002; Chen et al.2006; 
Delp et al.2009). It is an important differential 
expression technology and can be widely used 
to understand gene expression characteristics 
of organisms. Recently, DDRT-PCR has been 
successfully used to isolate expressed se- 
quence tags in a variety of organs. In this 
study, candidate genes of resistance to S. 
avenae were reported for the first time by us- 
ing this method. A new method of studying 
genetic resistance breeding by aiming at the 
same gene in all samples was proposed. 



Wang et al. 

Differential expression pattern of candidate 
genes of resistance to aphids 

Delp et al. (2009) reported that protease inhib- 
itors serine/threonine kinase exist in 
constitutive resistance of soybean cultivars 
resistant to Rhopalosiphum padi. This result 
implies that there should be some resistant 
genes in the resistant soybean cultivars such 
as 98-10-35. In this study, five kinds of differ- 
entially expressed patterns were described in 
the parents, F3 population, and backcross pop- 
ulation, which is a slightly different result 
compared to previous studies on cotton, 
wheat, and barley (Xing et al 2005; Meng et al 
2005; Zhang et al 2008; Zhu et al 2009; Zhang 
et al 2010). Four differential expression mod- 
els were detected between the hybrid and its 
parents in cotton: (1) Up expression only 
showed in hybrid and not in either parent; (2) 
Dominant expression showed in one of the 
parents but not in Fi and another parent, in- 
cluding the expression pattern in the female 
parent and the hybrid but not in the male par- 
ent, and the expression pattern in the male 
parent and the hybrid but not in the female 
parent; (3) The gene was silenced in one of 
the parents, including the expression pattern in 
the male parent but not in hybrid and female 
parent, and the expression pattern in the fe- 
male parent but not in the hybrid and the male 
parent; (4) Down expression showed in both 
parents but not in the Fi generation (Zhu et al 
2009). Previous studies also showed four dif- 
ferential expression models. In the present 
study, five expression patterns were found in 
the expression of resistant genes to S. avenae. 
At the same time, different bands in gene ex- 
pression were shown between parents and F3 
resistant populations and between the sensi- 
tive and resistant plants. It is best to select the 
target genes among lots of candidate genes by 
choosing from three patterns (2, 3, and 4 type) 
in the present study in order to save labor 
force and financial resources and overcome 
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the blindness of selection. This study is of 
significance for research on resistance of 
wheat to S. avenae and suggests a new ap- 
proach to study genetic resistance breeding. 
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Submergence Tolerant Using Differential 
Display Technique in Rice. Journal of 
Agricultural Biotechnology 14(6): 894-898 
(in Chinese with English abstract). 
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